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Anomalous negative heat capacities have been claimed as 
indicators of first order phase transitions in finite systems in 
general, and for nuclear systems in particular. A thermo- 
dynamic approach allowing for all Q value terms is used to 
evaluate heat capacities in finite van der Waals fluids and fl- 
nite lattice systems in the coexistence region. Fictitious large 
effects and negative heat capacities are observed in lattice 
systems when periodic boundary conditions are introduced. 
Small anomalous effects are predicted for small drops and for 
flnite lattice systems. A straightforward application of the 
analysis to nuclei shows that negative heat capacities cannot 
be observed for A > 60. 



The quest for discovery of the liquid to vapor phase 
transition in nuclei has progressed along two lines. On 
one hand, the cluster abundance has been studied as a 
function of mass and temperature and the phase dia- 
gram has been generated |^] in terms of Fisher's theory 
of clusterization On the other, caloric curves have 
been analyzed to look for plateaus and breaks possibly 
associated with the phase transition ; and fluctuations 
have been translated into heat capacities in the hope of 
discovering negative values which are widely considered 
to be indicators of phase coexistence in small systems 
i|. 

Regarding experimental caloric curves, it has been 
pointed out that their interpretation hinges upon an un- 
known pressure- volume relationship in the experiment, 
e.g. constant pressure or constant volume 

Regarding the heat capacities, negative values are typ- 
ically obtained in lattice gas, Ising, or Potts model cal- 
culations for finite systems ||8|"[l0|]. The origin of these 
negative values is generically attributed to the genera- 
tion of surface with its attendant energy cost, which is 
significant in small systems and which disappears in the 
thermodynamic limit. However, with numerical calcula- 
tions, it is not clear how this effect actually comes about, 
and if and how it may apply to actual systems like nuclei. 
The issue is all the more interesting since claims of nega- 
tive heat capacities have been made for nuclear systems 
i|. 

In this paper we are going to investigate the subject 
of caloric curves and heat capacities of finite systems in 
the coexistence region and the underlying role of varying 
potential energies ( "ground states" ) with system size on 
the basis of simple and very general thermodynamical 
concepts, in the hope of obtaining solid and unambiguous 



conclusions on these matters. 

Our study applies to leptodermous (thin skinned) van 
der Waals-like fluids and to models such as Ising, Potts, 
and lattice gas, which are capable of reproducing their 
general features. We shall show that: 

1. drops of leptodermous systems can indeed show 
negative heat capacities as a slight effect if allowed 
to evaporate at constant pressure; 

2. finite lattice systems can also present negative heat 
capacities as a slight effect if open boundary con- 
ditions are applied; 

3. finite lattices systems present negative heat capac- 
ities as a much greater effect if periodic boundary 
conditions are applied, the main signal being asso- 
ciated with the latter, very artificial conditions; 

4. nuclei cannot show negative heat capacities above 
A fa 60, while negative heat capacities are possible 
for A values below 60. 

Let us consider a macroscopic drop of a van der Waals 
fiuid with A constituents in equilibrium with its vapor. 
The vapor pressure p at temperature T is given by the 
Clapeyron equation 



AH„ 



dp 



(1) 



where AHm is the molar vaporization enthalpy and AVm 
is the molar change in volume. The Clapeyron equation 
gives a direct connection between what we might call the 
"ground state" properties of the system and the satura- 
tion pressure along the coexistence line. In fact, we can 
write 



AH,n = AE,n + PAV„, - AE„ 



T 



(2) 



and for T ^ AE, AHm ~ clv can be identified approxi- 
mately with the liquid-drop volume coefficient in the ab- 
sence of other terms like surface and Coulomb, etc. As- 
suming AVm ~ ym^°^ ^ that the vapor is an ideal gas and 
that AHm is constant with temperature, Eq. (|l|) can be 
integrated to give 



p ~ _Po exp 



AH„ 
T 



(3) 



Equation (|^) represents the p-T univariant line in the 
phase diagram of the system if AHm is assigned its bulk 



1 



value AiT^ . It was observed long ago that for a drop of 
finite size AHm must be corrected for the surface energy 



of the drop |11 

AHra = AH^^ 
= AHt 



Ai7f„ = AHl, - jSrn 

42/3 



(4) 



where 7 is the surface tension, Sm is the molar surface of 
the drop of radius r, as is the surface energy coefficient 
and is a geometrical constant. Substitution in Eq. (||) 
leads to 



p^Po cxp 



Pbuik exp 



AH" 

ri 

T 
a. 



(5) 
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At constant temperature the vapor pressure increases 
with decreasing size of the drop. In other words, each 
drop of constant radius r has its own r-dependent coex- 
istence line with the vapor as shown in Fig. 
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FIG. 1. The natural log of the saturated vapor pressure as a 
function of the inverse temperature for different droplet radii. 
The size of the open circles is proportional to the droplet 
radius. Arrows illustrate the path of evaporation at constant 
pressure. 

Let us now consider the case of isobaric evaporation 
of a drop starting from a drop with Aq constituents and 
evaporating into a drop with A < Aq constituents. Let 
us now define the drop size parameter 



y 



At constant pressure 



Ag-A 

An ■ 



(6) 



poexp 



T 



= po exp 



AH„,[yy 



Remembering that in most liquids |a^, 
in the lattice gas \ay \ = \as\, we have 



(7) 



and that 



AH^{y) 



from which follows 
Ty AH„,{y) 



1 



Al'\l-yYl\ 



Ai7£ 
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1 - 



Al'\\~yY/^ 



(8) 



(9) 



Thus, a slight decrease in temperature is predicted as the 
drop evaporates isobarically, thus leading to a negative 
isobaric heat capacity in the coexistence region as illus- 
trated in Fig. ||. In this figure, the abscissa is trivially 
related to Ai7, the heat absorbed by the evaporating 
drop at constant pressure, by the relationship 



AH = Aof dyAH„,{y) ^ AoyAH„ 
Jo 



(10) 



The same decrease can be more visually appreciated from 
Fig. (|^). As the drop is evaporating at constant pressure, 
the drop moves from one coexistence curve to another ac- 
cording to its decrease in radius, and thus to progressively 
lower temperatures. This slight effect is due not to an 
increase in surface as the drop evaporates, since the drop 
surface of course diminishes as A"^^^ , but to the slight in- 
crease of molar surface which does increase as A~^/^ as 
shown in Fig. ^. Also, the formation of bubbles in the 
body of the drop is thermodynamically disfavored by the 
factor / = exp(— 7AS'/T) where AS* is the surface of the 
bubble. 

It is worth pointing out that the smallness of the effect 
(~ 4% decrease in T as a droplet of A ~ 200 evapo- 
rates all the way to ^ ~ 100), is needlessly magnified 
by the widely practiced artful translation of the caloric 
curve into heat capacity, which jumps from -l-oo to —00 
as the slope of the caloric curve changes from constant 
to a slightly negative value. 

It might be argued that cluster formation may be re- 
sponsible for negative heat capacities. In physical fluids 
as well as in the Ising model the "mean" molecular weight 
of the equilibrium vapor remains very close to that of 
the monomer, as confirmed by the adherence of the va- 
por pressure vs. temperature to the Clapeyron-Clausius 
formula with constant AH„i Furthermore, if such 

an effect existed it would operate already for the infinite 
systems, which, patently, it does not. 

Let us now move to the amply studied cases of lattice 
gas, Ising, and Potts models. Here as above, the study 
of the evolution of the ground state properties can be di- 
rectly translated into the liquid-vapor coexistence prop- 
erties. We consider first an evaporating finite system in 
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FIG. 2. The temperature as a function of droplet size for a 
drop evaporating at constant pressure in a system with open 
boundary conditions. The sohd Una shows the case of a spher- 
ical drop, while the dotted line shows the case of a finite cubic 
lattice evolving as in Fig. ^ top. 



three dimensions of size = L^, with open boundary 
conditions. This case is essentially identical to the case 
of a drop discussed above (see Fig. ||). 

For maximal density at T = (the ground state) y — 
and the entire cubic lattice is filled. For decreasing 
densities, always at T = a single cluster of minimum 
surface is present, which evolves from a cube to a sphere. 
The associated change in surface is shown in Fig. (^). 

Since a cube has only a slightly larger surface than a 

sphere of the same volume by a factor of 6/ (^(Air)^^^ 3^^'^^ 
the resulting decrease of Ty at y = 1 is nearly exactly 
what a sphere of the same volume would experience in 
going from y = to y = 1/2 (Aq - 200). Thus the 
caloric curve from j/ = 0toy = l/2is essentially flat like 
in the infinite system, and the heat capacity is trivially 
infinite. 

So, where are the large effects reported in so many 
papers |p|~|Tot? We shall see below the effect created by 
the introduction of periodic boundary conditions. 

The introduction of periodic boundary conditions rids 
the system of "dangling bonds," as it were, by repeating a 
cubic lattice of side L periodically along the three coordi- 
nates. These conditions, originally introduced to mimic 
the infinite system, lead here to peculiar consequences. 

At y ~ 0, the lattice is filled with particles so that 
AH„i{0) = AiJ^j characteristic of the infinite system. 
As y increases at fixed lattice size, a bubble develops in 
the cube and surface is rapidly created. This is shown in 
Fig. (H). The bubble develops since the periodic bound- 
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y = (A^-A)/Ag 
FIG. 3. Top: The surface S° (dashed) and molar surface 
S!^ (solid) area of a drop for open boundary conditions nor- 
malized to their values at j/ = 0. Bottom: The surface 
(dashed) and molar surface (solid) area of a drop for peri- 
odic boundary conditions normalized to their values at y = 0. 
In-sets show the configurations at various values of y. 



ary conditions prevent 'as it were' evaporation from the 
surface. The bubble grows with increasing y until it 
touches the sides of the lattice. This occurs for y m 1/2. 
At nearly y = 1/2 and beyond, the "stable" configura- 
tion is a drop that eventually vanishes a,t y — 1. The 
change in surface associated with the range < j/ < 1 as 
well as the molar surface are shown in the bottom panel 
of Fig. I 

The evaporation enthalpy thus becomes 



from y — to y = 1/2, and 



y 



2/3 



A'J'il-y), 



AH^iy) ~ a„ 1 



Al^'{l-y) 



1/3 



(11) 



(12) 



from y = l/2toy = l. 

As a consequence, for periodic boundary conditions 



In 

Toe 



1 



y 



2/3 



<'(i-y) 



(13) 



from y = to y = 1/2, while from y = 1/2 to y = 1 
Eq. (|) holds. 

The results contained in equations (^ and (|l^) are 
very general and convenient. They bypass the clumsy nu- 
merical calculations for individual systems, and require, 
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FIG. 4. The temperature as a function of droplet size for a 
drop evaporating at constant pressure in a system with peri- 
odic boundary conditions. The sohd hne shows the case of a 
finite cubic lattice with periodic boundary conditions evolv- 
ing as in Fig. ^ bottom, while the dotted line and the dashed 
line are the same as in Fig. ^ and the vertical dash-dotted line 
indicates the case of 50% lattice occupation. 



for any case, just the knowledge of the geometry and of 
the surface energy coefficient. 

The dramatic effect of periodic boundary conditions 
can now be seen in Fig. ^ The temperature decreases 
substantially with increasing y, due to the fact that the 
molar enthalpy at ?/ = assumes its bulk value A_ff^ 
and must meet the previous case of open boundary con- 
ditions for y = 1/2. This may well explain the calcu- 
lated negative heat capacities reported in literature, as 
an artifact due to the unnatural choice of boundary con- 
ditions. The conclusion of this exercise is the following: 
if AHjniy) decreases with y (with decreasing drop size), 
for any reason, we expect negative heat capacities. Al- 
ternatively we expect the normal infinite heat capacities 
if AHm = constant or positive heat capacities if AHm 
increases with y. 

With the lessons learned above we can evaluate the 
heat capacities for nuclei. It is apparent from the above 
arguments that the key quantity is AHm and its depen- 
dence on the drop size, irrespective of the physical causes 
that determine its magnitude and dependence. In the 
case of nuclei the quantity AH^ is determined not only 
by the volume and surface, but also by all the other terms 
in the liquid drop model, such as the Coulomb and sym- 
metry energy all of which contribute to the mean binding 
energy per nucleon. Consequently one can immediately 
infer that when the binding energy per nucleon decreases 
with A, the heat capacity should be positive, and vice- 



versa. Thus, since the maximum binding energy per nu- 
cleon occurs at A ~ 60, negative heat capacities should 
be possible only for A < 60. Let us proceed more pre- 
cisely. We can rely again on the Clapeyron equations to 
calculate the heat capacity as follows 



Cp — 



but 



dH 


dH 


dA 




p^'dA 






dT 




dp 




dl 


p 


' dl 



dT 
dp 



(14) 



(15) 



From the integrated form of the Clapeyron equation we 
have 



dp 
dl 



1 dAH„ 
T dA 



so 



dT 
dA 

Finally 



1 dAHm TVm 
T dA ^ AHr, 



T dAH„ 

AHrr,. dA 



T 



~ ' dAH, 



(16) 



(17) 



(18) 



dA 



The derivative in the denominator can be evaluated ap- 
proximately from the dependence on the binding energy 
per nucleon B upon the mass number 



dAH„ 
dA 



dB 

Ia' 



(19) 



The liquid drop model allows us to estimate such a 
derivative. Without the Coulomb term, of course, we re- 
cover the results presented above for a drop: the binding 
energy increases with increasing A and tends asymptoti- 
cally to the value ~ 15MeV. Thus negative heat ca- 
pacities should be expected and the caloric curves should 
look like that shown in Fig. ^. The Coulomb and symme- 
try terms, however, become very important at large val- 
ues of A, say, along the line of /3-stability. From Fig. || it 
is apparent that the binding energy decreases with A for 
A >~ 60. Consequently in all this region of A, positive 
specific heats should be expected. Only for A <~ 60, 
negative specific heats are predicted. 

Despite the assumption of a temperature independent 
AH„i we expect these results to hold over a very broad 
range of temperatures. On one hand a monotonic change 
of AH,n with T should not alter the result; on the other 
hand it is known that for van der Waals fluids their vapor 
pressure p = p{T) is well defined up to the critical point 
by a constant AHm |l^ • Furthermore the appearance of 
clusters has minimal influence on AHm and therefore on 
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FIG. 5. Top: The binding energy per nucleon of atomic 
nuclei. Bottom: The associated heat capacity. 
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our results, because the vapor concentration is dominated 
by monomers ||^. 

This straightforward resuh based on elementary ther- 
modynamics and ground state binding energies raises se- 
rious questions as to the meaning of the negative heat 
capacities that have been claimed for large nuclear sys- 
tems. 

In conclusion, we have shown that: 

1. it is possible to generate caloric curves and heat ca- 
pacities in the coexistence region from the knowl- 
edge of the molar heat of vaporization, which must 
include all Q-value terms; 

2. simple drops and Ising, Potts models with 
open boundary conditions show minimal anoma- 
lies, while periodic boundary conditions introduce 
strong anomalies as artifacts; 

3. nuclei with A >~ 60 should present no anomalous 
negative heat capacities while this is possible for 
A <~ 60. 
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